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Abstract 

We demonstrate the common bihamiltonian nature of several integrable systems. The 
first one is an elliptic rotator that is an integrable Euler- Arnold top on the complex group 
GL(A, C) for any N, whose inertia ellipsiod is related to a choice of an elliptic curve. Its 
bihamiltonian structure is provided by the compatible linear and quadratic Poisson brackets, 
both of which are governed by the Belavin-Drinfeld classical elliptic r-matrix. We also 
generalize this bihamiltonian construction of integrable Euler- Arnold tops to several infinite- 
dimensional groups, appearing as certain large N limits of GL(A, C). These are the group of 
a non-commutative torus (NCT) and the group of symplectomorphisms SDif f{T^) of the 
two-dimensional torus. The elliptic rotator on symplectomorphisms gives an elliptic version 
of an ideal 2D hydrodynamics, which turns out to be an integrable system. In particular, we 
define the quadratic Poisson algebra on the space of Hamiltonians on depending on two 
irrational numbers. In conclusion, we quantize the infinite-dimensional quadratic Poisson 
algebra in a fashion similar to the corresponding finite-dimensional case. 
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1 Introduction 

One of the most interesting problems in the theory of integrable systems is a description of 
infinite-dimensional integrable systems with two or more space variables. In this paper we 
consider integrable hierarchies in (2-1-1) dimension, which we call elliptic and modified hydrody- 
namics. These are close relatives of the ideal hydrodynamics on a two-dimensional torus, which 
is known to be a highly non-integrable system. For instance, even the 4-vortex approximation 
of the ideal fluid dynamics is non-integrable [1]. We obtain these integrable versions of hydrody- 
namics starting from finite-dimensional integrable systems whose dimension of the phase space 
goes to infinity. Our approach is somewhat similar to the derivation of Toda field theory from 
the open Toda chain. In contrast with the Toda theory, however, we obtain non-local systems. 

The classical Euler (or, rather, Helmholtz) equation for the motion of an ideal fluid on the 
standard two-torus is 

dtS = {S,A-'S}, 

where S is the vorticity function of the fluid flow, { } is the Poisson bracket, and A is the Laplace 
operator on T^. The modified and elliptic hydrodynamics are defined on an elliptic curve, i.e., 
on a torus with a complex structure fixed. The corresponding equations are, respectively, 

dtS = {s,d-^s} 

and 

dtS = {s,p{d)s}, 

where d is the corresponding operator of the complex structure, and p is the Weierstrass p- 
function. (We postpone the precise description of p{d) till Section 5.) Note that to define 
the Laplace operator one needs to choose a metric, while for the operator 5 in the modified 
hydrodynamics is defined by a complex structure on T^. 

Such a modification of the fluid inertia operator from the Laplace operator A = dd, which 
depends on a metric on T^, to 9^ or p{d)^^, both of which depend on a complex structure 
on r^, brings in the intcgrability and even the bihamiltonian structure for the systems. We 
construct an infinite set of involutive integrals of motion with respect to two Poisson brackets. 
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One of the brackets is the standard Uncar Lie-Poisson brackets on the dual space to the algebra 
of the divergence- free vector fields on T^. The other Poisson structure is a quadratic Poisson 
algebra on Hamiltonians of vector fields. These two brackets are compatible and governed by 
the same classical r-matrix. Furthermore, we describe a recursion procedure for constructing 
the sequence of Hamiltonians for this linear-quadratic bihamiltonian structure of the hierarchy 
of the elliptic hydrodynamics, which thereby exhibits "the strongest form of integr ability." 

We come to this construction through the non-commutative deformation of to the non- 
commutative torus (NOT). The non-commutative deformation of the Lie algebra of vector fields 
is isomorphic to the Lie algebra of NOT. The latter is a special large N limit of gl(A^, C). We start 
with the elliptic rotators on GL(A^, C) [2, 3] and develop their bihamiltonian structure based on 
the Belavin-Drinfeld classical r-matrix [4]. The corresponding quadratic Poisson algebra is the 
classical limit of the Feigin-Odesski algebra [5]. We describe natural extensions of the elliptic 
rotators to the infinite-dimensional groups, which preserve their main properties. This brings 
us to bihamiltonian systems and integrable hierarchies on the NOT, which are of interest by 
themselves. 

As a byproduct, we also describe a quantization of the quadratic Poisson algebra on NCT. In 
this way we obtain an infinite-dimensional associative algebra with quadratic relations depending 
on the complex structure on and the Planck constant h &T^. 

Finally, the commutative (" dispersionless" ) limit of the elliptic rotators on NCT leads to the 
desirable bihamiltonian hierarchies of the elliptic and modified hydrodynamics. 

2 Main results 

1. The general setup [6, 7] 

Let G be a Lie group and g its Lie algebra. Consider an invertible linear operator J that maps 
the coalgebra g* to q. Its inverse operator J^^ is called the inertia tensor. The Euler-Arnold 
top corresponding to the group G is the Hamiltonian system on g* with respect to the linear 
Lie-Poisson brackets on g* and the Hamiltonian function given by the quadratic form 

i7 = -^(S,J(S)), Seg*, 

where ( , ) stands for the pairing between g and g*. Namely, the variation V-H" can be regarded 
as an element of g and the corresponding Hamiltonian equation of motion is as follows: 

dtS = {H,S} :=ad^^S. 

Recall that the Lie-Poisson brackets arc degenerate on g* and their symplectic leaves are coad- 
joint orbits of G. To descend to a particular coadjoint orbit O one should fix the values of 
Casimirs for the linear bracket. 

For some special choices of J the system becomes completely integrable. Some of the exam- 
ples are Manakov's tops on SO{N) [8], their limit N ^ oo found by Ward [9], the Korteweg-de 
Vries equation [10, 11], as well as the Camassa-Holm and Hunter-Saxton equations on the 
Virasoro group [12]. (In the infinite-dimensional case, the invertibility of J is understood as 
KerJ = 0.) 

2. EUiptic Rotators on GL{N, C) 

The elliptic rotator (ER) on GL{N, C) is an integrable Euler-Arnold top on this group, whose 
inertia operator is constructed with the help of the Weierstrass p-function of an auxilliary 
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elliptic curve, see Ref. [2, 3, 13]. Namely, let Tq,, a = {ai,a2) be the basis of sl(A'^, C) (B.4,) 
= (0, 1, . . . , A*" — 1), a 7^ (0, 0). The structure constants in this basis are 

Ce(a,/3) = ^sin7r0(a X e = k/N, 

where 1 < k < N and A:,A^ are coprime, see (B.8). For S = Y^a'^-aTa £ sl(iV,C)* the linear 
Poisson brackets assume the form 

{Sa,Sp} = Ce{a,(i)Sa+p. (2.1) 

Let p{x; r) be the Wcicrstrass function on the elliptic curve Ej- = C/(Z+tZ). Consider its values 
Pe(") = + Q!2r)0; r) on the lattice parametrized by a, the labels of the sl(A'", C)-basis. The 

inverse inertia operator is defined as 

3 : Sa ^ JaSa , Ja = P6»(a) , (2.2) 

It was proved in Ref .[2] that the equations of motion 

dtSo, = Y,So.-',Pe{l)S-y, (2.3) 

7 

defined by the Hamiltonian H = -i(S, J(S)) = 2'K'^ehi{^ ■ J(S)) and the brackets (2.1), have 
the Lax representation with the Lax operator U'°^{S,z) G s^A*", C) depending on the spectral 
parameter z G Et. The traces tr(L''°*(S, z))'^, k = 2,...,N being expanded in the basis of 
elliptic functions on E-,- produce the involutive integrals of motion and the hamiltonian H is 
among them. 

We call the commuting flows defined by these integrals the ER hierarchy. 

Let r{z), z e Et he the Belavin-Drinfeld classical elliptic r-matrix, see [13, 4] and Section 3 
below. The main result, which we will also generalize to the infinite-dimensional situation, can 
be formulated as follows. 

Theorem 2.1 i) In terms of the r-matrix, the linear brackets (2.1) can he written in the form 

{L{"\z),Ll°\w)}i = [r{z-w),Lr\z)+U2"\w)]. (2.4) 

ii) Consider the phase space extended by a new variable Sq, which has zero bracket with all the 
rest. The same r-matrix defines the quadratic Poisson algebra PN,e,T related to GL(A^, C) (3.19), 
(3.20): 

{Li{z),L2{w)]2 = [r{z-w),Lr{z)®L2{w)\, where L{z) = Sold + L''°\z) . (2.5) 

iii) The above two brackets are compatible, i.e. any linear combination of them is a Poisson 
bracket. 

iv) There exists a sequence of integrals of motion in involution with respect to each of the two 
brackets {hj,hk\i,2 = (here the lower indices refer to the linear and the quadratic brackets 
respectively). They provide the bihamiltonian structure of the elliptic rotator (ER) hierarchy 

{hj+i,S}2 = -{hj,S}i. 

The first two statements of Theorem for GL(A'^, C) are well known [5, 14]. Apparently, the 
bihamiltonian structure of the GLn elliptic rotators is new, and it gives the following 
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Corollary 2.1 The Casimirs with respect to one of the brackets generate non-trivial dynamics 
with respect to the other. 

In particular, the functional Sq, being the Casimir element of the linear brackets (2.1), leads to 
the equations 

dtS = {So,S}2, (2.6) 

that coincide with (2.3). 

It turns out that Theorem 2.1 also holds in the infinite-dimensional situation presented below. 

3. Elliptic Rotators on the Non-commutative Torus 

Consider the sin-algebra sing, called also the Lie algebra of the non-commutative torus (NCT). 
The Poisson brackets on the Lie coalgebra sm^ of the NCT has the form (2.1) with irrational 
number < ^ < 1 and the basis Ta = const ■ exp(27riQ;ixi) * exp{2TTia2X2) is parameterized by 
the infinite lattice a G Z©Z. Here exp(27ria:i), exp(27ria;2) are the generators of the NCT and * 
is the Moyal multiplication (see Appendix C). On the NCT we introduce the complex structure 
de^T depending on r (Im r > 0), and two real numbers e = {€±,€2), such that Oea are irrational 
and < Oca < 1 (see (4.3)). The inverse inertia operator J is the pseudo-differential operator 

3{s){x) = p{ed,,r)s{x) , 

where now S{x) = S-aTa{x). The equation of motion in the form of the Moyal brackets 
has the form 

dtS = {S, p{9d,,r)Sy , (2.7) 

where { }^ is the Lie-Poisson bracket on sing. In Proposition 3.1 we find the Lax form of 
(2.7). We construct the classical r-matrix (4.11) and prove the counterpart of the Theorem. The 
quadratic Poisson algebra P0^e,T (4.19), (4.20) on the NCT gives rise to the infinite bihamiltonian 
hierarchy, where the equations (2.7) correspond to the quadratic Hamiltonian functional with 
respect to the linear brackets. At the same time the equations (2.7) can be interpreted as the 
Hamiltonian equations (2.6) in the quadratic Poisson algebra P0,e,T with a linear Hamiltonian. 
For SL(2) this representation was found in [15] and for general case in [18]. 

4. Elliptic Rotators on SDiff{T'^) 

In the limit 9^0 the Lie algebra of the NCT becomes isomorphic to the Poisson algebra A of 
smooth functions on modulo constants Ham{T'^) ~ (7°°(T^)/C. The algebra of Hamiltonians 
can be also described by the corresponding Hamiltonian (or, divergence-free) vector fields on 
the torus. More precisely, the Lie algebra SVect{T'^) of divergence- free vector fields on T^ is 
the (universal) cocentral extension of Ham{T^), defined by the exact sequence 

^ Ham{T^) SVect{T^) ^ ^ , 

where the image of SVect{T'^) in is generated by the two fluxes eidi, €282. For i/j € HamiT'^) 
we have 

Vrii^) = ~^d2^ , V2{ii^) = ^diii^ , 

and V{^^^/j = \y{'il)),V{il))']. We construct the elliptic rotator on the Lie group SDif f{T'^) of 
area-preserving diffeomorphisms, corresponding to the Lie algebra SVect{T^). 
Consider the dual space of linear functionals Ham{T'^)* in the Fourier basis 

{e(a ■ x) := exp 27ri(a • x) \ (x = {xi,X2) , (a ■ x) = a\X\ + 02X2 , oij G Z} , 
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S = ^ Sae{-{a ■ x)) e Ham{T^)* . 

a 

The Poisson structure on A* assumes the form (cf. (2.1)) 

{Sa, Si3}i = (q X (3)Sa+l3 ■ 

This Poisson structure is degenerate and has an infinite set of Casimirs 

Ck= [ 8*=, where k = 2,3,... . 

On a coadjoint orbit O C Ham{T'^)* of SDiJfiT'^) the brackets are non-degenerate and values 
of the Casimirs are fixed. 

Define the operator J (after the rescahng) as the pseudo-differential operator p(9e^^), such 
that the Hamiltonian has the form 

H = ~J^^S-p{d,,r)S. (2.8) 

We prove that the corresponding Hamiltonian system admits an infinite set of commuting inte- 
grals, see Section 5. 

Recall that these elliptic rotators are parameterized, in particular, by the auxiliary elliptic 
curve. Consider the simplest version of the ER, corresponding to the rational degeneration of the 
elliptic curve Et- In the limit the (inverse) inertia operator becomes J = and the limiting 
Hamiltonian (2.8) takes the form 

see (5.23) and [15]. 

It is interesting to compare this with the Euler equation for an ideal fluid on a torus, where 
the inertia operator is the Laplacian (i.e., J = A"^ = (dd)^^), and the Hamiltonian is 

where S plays the role of the vorticity [7]. While the Euler hydrodynamics equation is highly 
non-integrable, the modification of the (inverse) inertia operator from Bd to 5"^ leads to an 
integrable hierarchy. We call the rational limit the modified hydrodynamics on T^. (In a sense, 
these systems are similar in spirit to the Etingof-Prenkel current algebras [16] on T^, which make 
use of the complex structure on the elliptic curve to construct a central extension of smooth 
currents on the torus.) 

Finally, we prove the Theorem 2.1 for the generic systems of the eUiptic hydrodynamics 
(2.8). In particular, the quadratic Poisson algebra 7^o,r (5.27), (5.28) on T^, along with the 
linear Poisson bracket, provide the bihamiltonian structure of the hierarchies of the elliptic and 
modified hydrodynamics. 

3 GL{N, C)-Elliptic Rotators 

3.1 Lax representation and integrals of motion 

1. The system description 

The elliptic GL(A^, C)-rotator is defined on the dual spaces of gl(A^, C) and sl(A^, C). Let S = 
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X^Q, S-^aTa, where Tq, is the basis of sl(A^, C), sec (B.4). Then the Poisson structure on the dual 
space sl(A^, C)* is given by the hnear Lie-Poisson brackets (B.8) 

{Sa, Sf3}l = C0{a, l3)Sa+/3 , (3.1) 
where Cg{a,P) are the structure constants of sl(iV, C) (B.8). The Hamiltonian has the form 

H'^ot ^ 27r202tr(S • J(S)) ^ -\Y1 S^p^S.^ . (3.2) 

7 

It defines the equations of motion 

atS = [J(S),S]. (3.3) 

Let J : Sa ^ JaSai be the (inverse) inertia operator with J a = pe{<^)-, where p0{a) are 
the values of the Weierstrass p-function and defined by (B.ll). Then (3.3) assumes the form 

dtSa = Sa-jS^P0i-f)C0{j, a) . (3.4) 



2. The linear brackets and r-matrix 



Proposition 3.1 (cf.[2]) The equations of motion (3.2) have the Lax form 

dtL'"\z) = [L'"\zlM{z)], (3.5) 

where 

L'ot = _^SaMa,z)Ta, (3.6) 

a 

M = -YSafe{a,z)Ta, (3.7) 

a 

and the functions ipe, fe are given by (B.13), (B.I4.). 
Proof. 

Substituting (3.6) and (3.7) in the Lax equation we obtain 

dtSa^e{a,z) = "^Sa-^S^iipeia - -f,z)fe{-f,z) - ipe{-f, z)fe{a - 7, z)) . 
7 

Now, using the explicit expressions (B.13) and (B.14) for ipg and fg respectively, as well as the 
Calogero functional equation (A. 22), we come to 

dtSa = Y Sa-ySjipg{a, z){po{a - 7) - peil)), 
7 

which coincides with (3.4). □ 

The Lie-Poisson brackets (3.1) admit the following r-matrix description [4, 13, 14, 17]. Define 
the classical r-matrix by 

rNfiA^ -w) = Yfe{l,z- w)T^®T_^ , (3.8) 
7 

where ips is defined by (B.13). 
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Lemma 3.1 [14]- The r-matrix (3.8) satisfies the classical Yang-Baxter equation 

+[^Ar,e,r(^), ^Ar,e,r(w^)] = . (3.9) 

Proof. 

We reduce the YB equation (3.9) to the following functional equation 

veil, z - w)ipe{a, z) - (pe{a + j,z- w)(pe{a, w) + (p0{a + 7, z)(p0{-^, w) = , 

by using the commutation relations in sl{N, C) (B.7). It easy to see that it can be rewritten for 
the functions (pg (B.12), since all three terms have the common exponent eg ((72 + (X2)z — 721^). 
Now, it coincides with the Fay identity (A.21), where we put 

ui = {"/ + a)6 , U2 = -jd , zi = z , Z2 = w. 

□ 

Proposition 3.2 In terms of the Lax operator (3.7) the brackets (3.1) are equivalent to the 
following relation for the Lax operator 

{L\"\z),Ll"\w)}^ = [rr,,eA^ - L'^'iz) ^Id + Id® r"\w)] , (3.10) 

Li = L^ Id, L2 = Id® L . 

Proof. 

To prove (3.10) we rewrite it in the basis (8) T^: 

{S-a, S-/3}(p0(a, z)ip0{P, w) 

= S_a-pC0{a, (3) {ip0{-P, z - w)ip0{a + /3, z) - (^^(a, z - w)ip0{a + /3, w)) . 
The same Fay identity (A.21) imphes that 

ip0{a,z)Lp0{J3,w) = ip0{-P,z- w)(p0{a + P,z)- (p0{a,z- w)(p0{a + P,w) . 
Thus we come to (3.1). □ 

3. The hierarchy of the Lax equations 

The Lax operator L'^'"* (3-6) has the following properties: 

i) L^°^ is an sl(A^, C)-valued meromorphic function on Et with a simple pole at the origin 
satisfying 

Res L^"\z)U=o = J2S-aTa. 

a 

ii) L**"* satisfies the quasi-periodicity conditions 

L'-^^z + 1) = QL'"'\z)Q-^ , U°\z + t) = kU°\z)A-^ . 

These properties imply that Tr{L{z))^ are doubly periodic functions with the poles up to the 
order k. Thereby, they can be expanded in the basis of the Weierstrass function and its deriva- 
tives 

ir{U'^\z)f = Io,k + L2,kP{^) + ■■■ + Ik,kP^'"'\z) . (3.11) 
In particular, in this way we obtain the Hamiltonian 

tr(L-*(z))2=/o,2 + /2,2p(^), Io,2 = 2(-^) H, /2,2 = trS^ 
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Proposition 3.3 The coefficients Ig^k involution with respect to the linear bracket (3.1) 

{Is,k,Im,3}l=^- (3-12) 

Prooi 
We have 

ML\"\z)fMLl"\w)y}, = tr{(Lr*(z))^ {Li"\w)y}, . 

Then, it fohows from (3.10) that these functionals Poisson commute. Using the expansion (3.11) 
we arrive at the involutivity of the coefficients (3.12).n 

In particular, all functions Ig^k Poisson commute with the Hamiltonian H (3.2). Therefore, 
they play the role of conservation laws of elliptic rotator hierarchy on Gh{N,C). We have a 
tower of ^^^^^^ independent integrals of motion 

Io,3 l2,3 h,3 

Io,n h,N In,N 

Note that Ik^k: k = 0,2,3 ... ,N are the Casimirs corresponding to the coadjoint orbit 

7^-* = {SGgl(A^,C), s = 5-'s(%}. 

The conservation laws Ig^k generate commuting flows on 71'^°^ 

ds,kS = {Is,k, S}i , {ds,k '■= dt,^k) ' 

that wc call ERn hierarchy. 

In what follows we will need another set of the conservation laws coming from the coefficients 
of the spectral curve 

C : F{X,z)=det{X-Id + L'^''\z)) = 0. (3.13) 

Here 

F{\, z) = X^ + X^-'^a2{z) + ... + aN{z) , 

ak{z) = Jok + J2kp{z) + ... + JkkP^^~'^\z) . (3.14) 

Since ak{z) and Tr{L{z)Y are related by the Newton formula, one can recursively represent Jg^k 
as polynomials of Im,ji j k, m < j 

Js,k = Y.''tX] n {i = {hM,...)). (3.15) 

For example, 

Js,2 = -Is,2, {s = 0,2); Js,3 = -Is,3, (s = 0,2,3); 

Jo,4 = ^/o,4-/o'2-§^l,2, (52 = 60^(m + nr)-^); (3.16) 

m,n 

J2,4 = --^2,4 — 2/0,2-^2,2 ; J3,4 = 2^3,4 ■^4,4 = ^-^M " 1^,2 i 

where in the second line we have exploited the relation 12p^ = 2p" + 5(2. 
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3.2 The classical GL{N, C) Feigin-Odesski algebras 

Here we consider the Feigin-Odesski Poisson brackets related to GL(A'^, C). In notations of 
Ref. [5] it is the classical limit of the quadratic associative algebra Q^^^k- 

It turns out that the same r-matrix (3.8) defines a quadratic Poisson algebra. We start with 
the Lax operator (3.6). Modify it in the following way 

L{z) = -SoId + U'°\z), (3.17) 

where 5*0 commute with Sa with respect to the linear brackets (3.1). Define the brackets between 
the entries of L as follows 

{Li {z) , L2 (w) }2 = [rN,e,T (z-w), Li {z) ® L2 {w)] , (3.18) 

(see Ref. [13, 14]). These brackets are Poisson, since the Jacobi identity is provided by the 
classical Yang-Baxter equation (3.9). A finite-dimensional quadratic Poisson algebra is extracted 
from (3.18) by eliminating the dependence on the spectral parameters z,w. 

Proposition 3.4 The quadratic Poisson algebra on gl(iV, C)* has the form 

{Sa, So}2 = ^ Sa-jS^ipei-y) - pg{a - 7))Ce(a, 7) , (3.19) 

{Sa,Sp}2 = SoSa+pCe{a,(3)+ J2 Sa-jSp+^fg{a,(3,j)Ce{7,a-P), (3.20) 

where Cg{a,P) are the sl(iV, C) structure constants (B.8), 

fe(a, /3, 7) = Ceil) + CeiP - a + 7) - Ce(/3 + 7) + Ce{a - 7) 

and -de, Q are the zeta constants (B.IO), (B.ll). 
Proof. 

The equation (3.18) yields the following form of the brackets 
{Sa, Sb}ipe{a, z)ipe{b,w) = '^Sa-jSb+jipe{j,z-w)ipe{a-'y,z)ipe{b + -f,w)Ce{'y,a-b) . (3.22) 

Note first, that the common exponents in the expressions of the functions cpQ (B.13) coincide in 
the left and the right sides. This allows one to pass to the functions (pQ and in this way to use 
the Fay-type identities. We rewrite the last expression in the form 

{Sa, Sb}(t)e{z)(t)0{w) = ^ 5'a_^5'b+^C0(7, a-b) 

7 

X {4>e{l, z - w)(j)0{a - 7, z)(j)e{b -I- 7, w) - 4)eia - 6 - 7, z - w)(t)e{b + 7, z)(t)e{a - 7, w)) . 

Consider first the case 6 = 0. Then we can apply (A. 30) for u = (ai + a2T)9, v = (71 + 72)^- It 
immediately leads us to (3.19). If a, 6 / we have the terms of two types. For 7 7^ a, —6 one 
uses (A.27), (A.28) and finds the second term in the right side of (3.20). If 7 = a or 7— = b one 
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should take into account that (l)g(0,z) = 1 and then use the Fay identity (A. 21). This gives us 
the first term. The identity (3.21) follows from (A. 29). □ 

We denote by Pn.o.t the quadratic Poisson algebra (3.19), (3.20). Recall that r is the modular 
parameter of the auxiliary curve Er and 6 = k/N. The algebra - r classical Sklyanin 

algebra [14]. In this case the r.h.s. of (3.20) contains only the first term. 

(2) (2) (2) 

As for the linear brackets, there exist N Casimirs Cg ,6*2 , . . . , C}^ of the algebra PN,e,T 
[18]. They can be read off from the expansion of the spectral curve (3.13) in the basis of elliptic 
functions 

N 

det(So • Id + V°\z)) = C^') + Ck^P^^'''\^) ■ (3-23) 

fc=2 

It follows from (3.13), (3.14) and (3.17) that they are linear combinations of the integrals Jg^k 

N 

Cq^^ = Sq + ^ Sq "^Jo,m , (3.24) 

m=2 

JV-fc 

^ = E Sf"'Jk,k+m , ik = 2,...,N). (3.25) 

m=0 

3.3 The bihamiltonian structure 

Two Poisson structures are called compatible (or, form a Poisson pair) if their linear combina- 
tions are Poisson structures as well. 

Proposition 3.5 The linear and quadratic Poisson brackets on gl{N, C) are compatible. 
Proof 

Let us replace 5*0 Sq + X in the quadratic brackets (3.19), (3.20). Then (3.19) does not change, 
while (3.20) acquires an additional linear term. In this way we define a one-parameter family of 
the brackets 

{S,S}a:={S,S}2 + A{S,S}i, (3.26) 

where {S, S}i denote the linear brackets (3.1). Therefore, the linear combination {S, S}a of two 
brackets can be obtained from the quadratic bracket by a simple shift of Sq. The result of the 
shift, of course, still satisfies the Jacobi identity, and hence it is a Poisson bracket for any A. □ 

We denote this family of the quadratic Poisson algebras by PN,e,T,\- The algebras are iso- 
morphic for different A and degenerate to the linear Poisson algebra gl{N, C)* as A ^ oo. 
Consider the Casimir functions of the Poisson algebra PN,e,T,x 

%(A) = cf (So + A) {k = 0,2,...,N). 
It follows from (3.24) and (3.25) that /ife(A) are polynomials in A 

N-k 

hk{X) = Y hm,k\"' . (3.27) 

m=0 
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The coefficients are defined as 
It implies that 



, . A=o 

a! 



h - C^^-" ^ \^ {N-m)\ N-m-a J .o HON 

{N — a)ia\ ^ [N — m — ay.ai 



m=2 

N-k-a 



^».^= E ^ go^-^-"^-Vfc,feH-r. , (fe = 2, ■ ■ ■ , TV) . (3.29) 

In particular, 



/i^_i,o = iV5o , = ^^^^ ^^ ^g + -/o,2- (3.30) 

Conversely, one can express Jm,k+m as a linear combination of hN-m-s,m-i s = 0, . . . ,k from the 
relations (3.28) and (3.29). For example, 

Jjj = hN~j,j , Jj,j+i = hN-j-ij -{N - j)SohN-j,j , 

T h (AT -^Qh (iV_j-)(iV_j_l) 
Jj,i+2 = tlN-j-2,j - (iV - J j^o/in-i-l,i H ^ b^tlN-jJ ■ 

We arrange the quantities k in the triangular tableau 



^0,0 ^1,0 

^0,2 ^1,2 

^0,3 ^1,3 



hN-3,2 hN-2,2 
hN-iA hN-3,3 



0,N-1 



h 



ho. 



l,N-l 



N 



(3.31) 



Note, that the second line corresponding to /i^ i is absent since hi{\) = 0. The left side 

(2) 

of the triangle contains the Casimirs hoj = C- (3.24), (3.25) of PN,d,Ti while the right side 
represents the Casimirs hi\f-jj = Jjj of the linear brackets on gl(A^, C)*. 

The remarkable property of the quantities ha^^ is ttiat they are in involution with respect to 
both the linear and quadratic brackets: 

{ha,k, hb,j}l,2 = , 

see, e.g., [19]. This fact follows from the identities 

{hk{X),hj{\)}x = 0, (3.32) 

held for each A. 

Proposition 3.6 The integrals ha^k provide the recurrence relation for the elliptic rotator hier- 
archy ERn 

{ha,k, S}i = — {/ia+l,fe) S}2 . (3.33) 
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Proof. 

Since /ifc(A) are the Casimirs, we have {/ifc(A), S}\ = 0. Substitute the representation (3.27) into 
this equation. Then the recurrence relation (3.33) comes from the coefficient in front of A". □ 

In this way one can start with the Casimirs of the hnear brackets and produce a non-trivial 
dynamical system using the quadratic brackets 

{Jk,k, S}2 = -{Jk,k+i, S}i , (3.34) 

(see (3.30)). In particular, the flow (3.4) corresponding to H = 7o,2 can be represented by the 
quadratic brackets with Sq: 

dtS = {So,s}2. 

Alternatively, one can start with the Hamiltonians in the left side of the table (3.31) and the 
linear brackets. 

This allows us to conclude the proof of Theorem 2.1 for the hierarchy ERjv. Namely, Propo- 
sitions 3.2, 3.4, 3.5, and 3.6 are equivalent, respectively, to i), ii), iii), and iv). 



4 Elliptic Rotators on non-commutative torus 

4.1 The Lax equation and integrals of motion 

1. The system description 

The elliptic rotator on the non-commutative torus As is a generalization of the elliptic rotator 
on GL(A'", C). We consider the Euler- Arnold top on the group SINg of the NCT algebra A$ 
(Appendix C). The Lie coalgebra sirig is equipped with the linear Poisson brackets 

{Sa,Sp}i = Ceia,P)S^+p, (aGZ(2)), (4.1) 

where S = ^ S^aTa and 

Z(2) = {a = (ai,a2), aj e Z , 0^(0,0)}. 

The inverse inertia tensor J maps sirig sing, and depends on four parameters: 6,t G 
C, {Quit > 0) and e = (£1,62)) where < < 1 and e^^ are irrational numbers. The 
components of J are given by the elliptic constants Jq, = pe,e(o) (C.18) 

J(S)= P0M)S-aTa- 

The Hamiltonian has the form 

H = 27r^9^ [ S-J(S) = -J ^ p,,,(a)5a5_,. (4.2) 

Introduce a complex structure on the NCT Ae depending on e = (ei, €2) and r such that for 
X = J2a (^aTa One has: 

Be^rX = ^(eioi -I- e2a2r) (4.3) 
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Then, in the Moyal representation, the operator J can be identified with the pseudo-differential 
operator 

3{S){x) = p{ed,,r)S{x), (4.4) 

and 



f S-p{ed,,r)S. (4.5) 



'Ag 

The equation of motion in the form of the Moyal brackets has the standard form 

dtS = {S, ip{ed,,r)S)}' := ad*(s)S , (4.6) 

or, in the components, 

dtSa= S^pe,e{l)Sa-j . (4.7) 

7ez(2) 

2. The r-matrix and the Lax equation 
Define formally the Lax operator 

L'-"\z) = - J2 S^ipeAa,z)Ta, (4.8) 
aez(2) 

where ipe^e is given by (C.19). Below we formulate the conditions on the phase space, which 
show that L^°^{z) is well defined for z G -Er, z ^0. Note that 



Ja„ 



0. 

lAg 

The following Proposition is an infinite-dimensional analog of Proposition 3.1. 
Proposition 4.1 The equations of motion (4-7) have the Lax form 

dtU°^ = [L™*, M™*] := {U'''\ M™*}^ (4.9) 

with 

M'"\z) = - S^fe,e{a^z)T^, (4.10) 

Qez(2) 

and fg^^ (C.21). 

The proof is analogous to the finite-dimensional case. It is based on the Calogero functional 
equation (A. 22). 

Introduce the classical r-matrix on sing ® sing as the following sum: 

rg^eA^ - = Y feAl^Z -'^)T'y®T_-y, (4.11) 

7ez(2) 

Lemma 3.1 and Proposition 3.2 have the following analogs for the NCT: 

Proposition 4.2 a) The r-matrix (4-11) satisfies the classical Yang-Baxter equation (3.9). 
h) The r-matrix defines the Poisson brackets for the entries of the Lax operator (4-8) 

{Ll°\z), Ll°\w)}i = [rg^eA^ - w), L™*(z) (g) /d + /d (g) L^''\w)] . (4.12) 

They are equivalent to the linear brackets (4-1) on sin^. 
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The proof of Proposition 4.2(a,6) is similar to those of Lemma 3.1 and Proposition 3.2 and is 
based on the Fay identity (A.21). 

3. The hierarchy of the Lax equations 

In order to construct integrals of motion we first discuss the properties of the Lax operator. 
This operator is the meromorphic quasi-periodic function on E-,- with a simple pole 

i?esL^°*(z)U=o = S, 

L™*(z + 1) = C/i^iL™*(z)C/f ^1 , + r) = U^' L'-''\z)U^'^ . 

where Ui, U2 are the generators of NCT Ae (C.l). (This follows from (4.8), (A.9) and (B.15).) 
It means that the integrals 

are doubly periodic functions on Et with poles up to order j. Thereby, we have the expansion 




(4.13) 



In particular. 



Note that 




are the Casimirs of the linear brackets. 

The phase space TZ of the elliptic rotator is determined by the following properties of the 
integrals Isj 

n = {SeA*g\i) Is j < 00, ii) lim < 1 for all s < 3} . (4.14) 

The first condition means that the traces (4.13) of the Lax operator are well defined on Et for 
z ^Q. In particular, we have 

/o,2<oo, i-c. ^ (a) S'aS'-Q, < 00 . (4-15) 

We will use the second condition below to define the Casimir functions for the family of quadratic 
Poisson brackets. 

Further, due to Lemma 3.1, an analog of Proposition 3.2 is the following. 
Proposition 4.3 The quantities Ig^k pairwise Poisson commute. 

The conservation laws Ig^k generate commuting flows on the phase space TZ with respect to the 
linear brackets 

ds,kS = {Is,k,S}'i , (4.16) 

where dg^k stands for the corresponding time derivative dt^ ^ ■ We call these equations ERq 
hierarchy on the NCT. 

One can show that all flows can be represented in the Lax form with L^°* (4-8) and the 
corresponding Mg^k- Furthermore, there exists a set of the integrals of Js^^-type. They are 
related to the integrals Imj by the same formulae (3.15), (3.16). 
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4.2 Quadratic Poisson algebras and the bihamiltonian structure on NCT 

To define the quadratic Poisson algebra on the phase space Tl (4.14) we modify the Lax operator 
(4.8): 

L = -5o • To + L™* = So-To + Y, S-a^eA^^^ z)T^ > (4-17) 

a 

where 5*0 commutes with Sa with respect to the linear brackets. Due to the YB equation, the 
r-matrix (4.11) defines the quadratic Poisson brackets 

{Li(z), L2{w)}l = [rg^eA^ - w), L^{z) ® L^iw)] . (4.18) 

Again, we can get rid of the dependence on the spectral parameter. 

Proposition 4.4 The quadratic Poisson brackets on the phase space TZ are defined as follows 
{Sa, S0Y2 = J2 Sa-jS^PeAl) - PeA^ - 7))Ce(«,7) , (4.19) 

{.5a,5;3}^ = 5o5«+;3Ce(a,/3)+ .5a-75/3+7fe,e(«,/3,7)Ce(7,«-/3), (4.20) 

77^a,-/3 

where 051(0,7) is a structure constant of the sing-algebra, 

fe,6(a,/3,7) = Ce,e(7) + CeAP - « + 7) - CeAf^ + 7) + CeA» " 7) (4-21) 

^^(0)^e,,(a)T?e,,(/3)^e,,(/3 - a + 27) 
i?e,.(a - l)^eAP + l)^9AP " " + l)^9Al) ' 
and 'do are the constants (C.17), (C.18). 

Proof 

The equality (4.18), being reduced to the coefficients in the front of Ta ® T^, can be rewritten 
in the form (3.22) where the functions ifeia, z) are replaced on (feA'^^ Then using the cubic 
functional equations (A. 27) and (A. 28) wc obtain to the algebra (4.19), (4.20). Due the definition 
of the phase space (4.15) the series in the right hand sides of (4.19), (4.20) converge. □ 

We denote this Poisson algebra by Pe,e,T- Set now ei = £2 = 1. The corresponding algebra 

^e,i,T can be considered as a special large N limit of the finite-dimensional algebras PN,e,T (3.19), 
(3.20). For this we just replace the rational number 6 = k/N by an arbitrary irrational number 
< ^ < 1 in the algebra Pn,0,t- 



Proposition 4.5 a) The linear (4-1) o,nd quadratic (4-19), (4-20) Poisson brackets defined on 
TZ are compatible. 

b) The ERs hierarchy admits the bihamiltonian structure. 
Proof. 

One can shift + A and define the family Poisson algebras Pe,e,T,\- To define the Casimir 
elements on PN,e,T,\ we consider 

logdet(So + A + L-*(z)) = log(5o + A) + .^—>.t,L{zf . 
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Then we come to the infinite set of Casimirs ^ 

/ 1 N fc+l 



k=2 



k=2 

These functionals are well-defined on the phase space (4.14). This allows us to introduce the 
new set of conserved quantities 



s 



The latter Poisson commute with respect to the both types of brackets and give rise to the 
bihamiltonian structure of the hierarchy 

{hsj,S}'i = -{hs+ij,S}l (4.22) 

In particular, we can represent the flow (4.6) in the form 

dtS = {So,S}l 

(see (3.30)). □ 

This proposition concludes the proof of Theorem 2.1 for the hierarchy ERg. 



5 Elliptic Rotators on SDiff{T'^) 
5.1 Description of the hierarchy 

In the " dispersionless" limit 6^0 the Lie algebra sing turns to the Lie algebra Ham{T^) of 
Hamiltonians on a two-dimensional torus, see (C.22). This algebra of Hamiltonian functions can 
be represented by the Lie algebra of the corresponding divergence-free vector fields SVect{T^). 
More precisely, to pass from Ham{T'^) to SVect{T'^) one has to discard the constant Hamiltoni- 
ans, but add the "flux" vector fields d/dx\ and d/dx2 corresponding to multivalued Hamiltonian 
functions xi and X2 on the torus. 

We define the elliptic rotator system ER on the Lie group SDiff{T'^) by considering the 
limit ^ — of the ERg-system described above. Let 9^0, e\^2 oo, such that 

lim(0ei,2) = e'i,2 < 1 ; ^1,2 irrational. (5.1) 

In what follows wc drop the superscript '. 

Let S = S^ae{a-x) G Ham*{T^), where e{a-x) is the Fourier basis (C.23) ofHam*{T^). 
In the Fourier basis, the linear Poisson bracket assumes the form 

{Sa,Sp}i = (a X P)Sa+p , axp = aih - Q;2/?i ■ (5.2) 
^Since we pass from det to logdet the basis of the Casimir functions differs from (3.24) 
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The inverse inertia operator J : Ham*{T^) — > Ham{T'^) becomes 

3 : pe{ot)Sa , pe{a) = p(eiai + e2a2T; r) , 

where a G Z*^^) (B.5). The operator is well defined since ej are irrational. In other words, the 
operator J is the pseudo-differential operator 

J : S(x) ^ p(4,r)S(x) , 

where Sg^r is the operator of the complex structure on the elliptic curve, commutative torus T^. 
In fact, the complex structure depends on the ratio re2/ei. 
The quadratic Hamiltonian of the system is 

^ = E SjPe{l)S-y = -\ J^^ S{pid,,r)S) . (5.3) 

and the corresponding equations of motion are 

dtS = {S,p{de,r)S}i. (5.4) 

This equation is an elliptic analog of the hydrodynamics equation on the torus , regarded as 
an elliptic curve, see Remark below. We call it the elliptic hydrodynamics. 
Define the Lax operator 

L''"\x;z) = - ^ Sa<Pe{a, z)e{a ■ x) , ip^{a,z) = ip{eiai + Te20i2,z) . (5.5) 

aGZ(2) 

The conditions on the phase space formulated below, see (5.12), ensure that the operator 
L^°*(x, z) is well defined for z G E^-, z 7^ 0. Note that 



/ i^r*(^) = o. 



Proposition 5.1 The equations of motion (5.4) have the dispersionless Lax representation 

with given by (5.5) and 

M'''\x;z) = - J2 Safeia,z)eia-x), (5.6) 

where 

fe(a, z) = /(eitti + Te2a2, z) = e{e2a2)du(j){u, z)\u=eiai+re2a2 ■ 

The proof of Proposition 5.1 is the same as Proposition 3.1. 

Note that operators L^"^ and M^°* satisfy the quasi-periodicity properties 

L'-°\xi,X2;z) = L^°\xi + e2,X2;z + l), (5.7) 

L'-°\xuX2; z) = L^''\xi,X2 - ei; z + t) , (5.8) 

M'°\xi,X2; z) = M^^^xi + e2,X2;z + l), (5.9) 
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and 

M™* {xi,X2;z)- M™* (xi , X2 - ei ; z + r) = 2^iL™* {xi,X2,z). (5.10) 
Furthermore, it follows from (5.7) and (5.8) that there exists the expansion 

k 

[ (L™*(x,z))'^ = /o,fc + ;^W^-2)(z). (5.11) 

The phase space 7iP of the elhptic rotator on SDiff{T^) can be described similarly to the 
phase space TZ for SINg-group, see (4.14): 

7^° = {S e Ham*{T^) | i) < oo , ii) lim < 1 for all s < j} . (5.12) 

Define the r-matrix on (8) by 

reA^^y^^) = XI </'e(a,2:)e(Q;-a;)e(-a-y), z G (5.13) 

Proposition 5.2 a) T/ie r-matrix (5.13) satisfies the "dispersionless" YB equation 

+{r^A^, v; z), r^^r{y, v, w)} = . (5.14) 
b) In terms of the Lax operator (5.5) the canonical brackets on take the form 

{L-*(x; z), L-*(y; w)} = K,(x, y, z - w), U"\x- z)} + {r(x, y, z - w^'^^y; w)} . (5.15) 

Proof 

These statements are consequences of the Fay identity (A.21), cf. Proposition 4.2. □ 
The form of the brackets (5.15) implies that 

{/ (L™*(x,z))^ / (L-*(x,z))n = 0. 

In turn, then (5.11) produces the infinite sequence of conservation laws Is,k in involution. They 
define the hierarchy of commuting flows 

ds,kS = {Is,k,S}i, (5.16) 

which is the dispersionless limit of the hierarchy (4.16). The equations can be represented in 
the form of the dispersionless Lax equations. 

5.2 Quadratic Poisson algebras, and the bihamiltonian structure on SDiff{T'^) 

To describe the quadratic Poisson brackets, we first pass from the Lie algebraic Lax operator 
L^°* to the operator 

L, = -So + L'°\x;z). (5.17) 
Define the quadratic Poisson algebra by the formula 

{Le{x; z), Le{y; w)} = {re^^, y;z- w),L^{x; z)L^{y; w)} . (5.18) 
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Proposition 5.3 In terms of the Fourier modes, the Poisson algebra (5.18) has the form 

{S'a, Sq}2 = ^ Sa-^S^{pe{l) - Pe{a - 7))(a X 7) , (5-19) 

{S^,Sp}2 = SoS^+f){axP)+ 5,-7'5/3+7f.(a,/?,7)(7X (a -/?)), (5.20) 

where 

fe(a,/3,7) = Ce(7) + Ce{f3 - a + 7) - C(/3 + 7) + Ce{a - 7) 
^ T?^(0)T?,(a)T?,(/?)z^,(/?-a + 27) 

Ceil) = C(ei7i + ^£272), Ail) = C(ei7i + ^€272)- 

The proof of this Proposition is similar to that of Proposition 3.2. 

By the same trick as before we can show that the hnear and the quadratic brackets are 
compatible. It allows us to construct the integral of the hg^k type and define the recurrence 
representation (4.22) for the bihamiltonian structure on the hierarchy. In particular, the equation 
(5.4) can be represented as 

dtS = {So,s}2. 

5.3 Rational limit of the elliptic hydrodynamics 

1. Description of the limit 

So far we have been dealing with the elliptic curve Et parameterized by the two half-periods 
i^i,i^2, T = LO2/UJ1. WE replace our main functions in the following way: 

Ce(7) = C(ei7i'^i + e272t^2;t^i,a;2,) , peij) = p(ei7i'^i + ^272^2; wi, a;2, ) . (5.21) 

Now we are going to consider the rational limit of the elliptic curves, i.e., degeneration lima;i^2 
00. We look at the double scaling limit 

limei^2^0, lima;i^2 00 , limeict'i = l, lime2u;2 = r. 

Then lima;i^2 00 leads to the rational degeneration of the elliptic functions (5.21) 

lim Ce(7) = — > lim Peil) = \ ^ 7r = 7i + 72^ 

wi,2— »oo aii,2— >oo -y^ 

This implies, in particular, that the inverse inertia tensor assumes the form 

J(S) = a-2S(a;i,X2)= V \S.^T^, (5.22) 

where = Q;i + C)i2T- The equations of motion are defined by the corresponding quadratic 
Hamiltonian 

H = -\ f Sr^S = - J V ^^a^-a . (5.23) 

with respect to the linear Poisson bracket on Ham{T'^): 

dtS = {S,d-^S)}i. (5.24) 
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In the Fourier components the latter becomes 



dtSa — ^ ^ S^—Sa—j ■ 
7ez(2) 

We call these equations the modified hydrodynamics on the torus . 



Remark 5.1 Recall that the standard hydrodynamics of an ideal fluid on the torus is given 
by the Euler equation 

dtS = {S,A-'S}i, 

on the vorticity function S. The definition of the Laplace operator uses metric, while to define 
the modified hydrodynamics (5.24) we need to fix a complex structure on . 

Now consider the Lax representation for the modified hydrodynamics. In the double scaling 
limit we have 

liraipJa, z) = ex.p(2'Kiaiz) I 1 — ) , 

lim fJa, z) = — exp(2Triaiz)—^ 
Therefore, the Lax pair (5.6) in the rational limit is 

L{xi,X2;z) = ~'^SaCxp{2TTiaiz) ( \- -jTa- 

M(xi,X2) = — 5a exp(27riQ;i2;)— ^tTq, . 

a ' 

We can drop the exponents in the both operators and come to the following expressions 

L{xi,X2;z) = —d~^S{xi,X2) — ^S{xi,X2) , (5.25) 

and 

M{xi,X2) = -d-^S{xi,X2) . (5.26) 
It is easy to see that the dispersionless Lax equation 

dtL{xi,X2;z) = {L{xi, X2; z), M{xi, X2)} 

is equivalent to the equations (5.24). 

The infinite set of the integrals of motion comes from the decomposition 



(-1)^ / L^{xi,X2]z)dxidx2 = Io,k + 



s=2 

The corresponding infinite hierarchy of the modified hydrodynamics 

ds,S = {Is,k, S}l 
has the dispersionless Lax form representation. 
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2. Bihamiltonian structure 
It turns out that the bihamiltonian structure survives in this hmit. In fact, we have the fohowing 
quadratic Poisson algebra Vo,t on the commutative torus 

{Sa, So}2 = Sa-^S^{^ - ^ ){a X 7) , (5.27) 

{Sa,Sf3}2 = SoSa+/3{ax(3)+ ^ Sa-jSp+^f{a,P,j){jx{a-(3)), (5.28) 

where 

/(a,/3,7) = ( — + ^— — ^ — + ), 7^ = 7i+T72. 

Similarly to the above, we can consider a one-parametric family of the quadratic algebras by 
replacing Sq ^ Sq + A. Note that it degenerates to the standard Poisson brackets on in 
the limit A — > 00. This allows us to define the bihamiltonian structure for the hierarchy of the 
modified hydrodynamics. In particular, the equation (5.24) can be rewritten in the form as 

dtS = {So,s}2. 

Such a bihamiltonian structure is a curious feature in the modified hydrodynamics, emphasizing 
its drastic difference from the classical hydrodynamics. 



6 Quantum counterparts 

In this section we present two associative algebras "quantizing" the discussed above two Poisson 
algebras in the case of the NCT. It is easy with the linear bracket. Indeed, the quantization of 
the Lie- Poisson algebra on sirig leads to the universal enveloping algebra sing. 

Replace the classical variables S on the NCT Ae by non-commuting variables S. Now we 
consider the quantization of the quadratic Poisson algebra Pd^e,T ^-nd construct the associative 
algebra Pg,c^T.hi where fi G Er is the deformation parameter. For this wc introduce the quantum 
i?-matrix related to the group SINg. The quantum Yang-Baxter equation defined on GL(A'^, C) 
[13, 14] is generalized to SINg in the following way. The quantum i?-matrix on SINg SINg 
assumes the form 

R{z, w) ¥'e,e(c, z - w\h)Tc ® r_c , 

c 

where Tc is the basis of SINg and the main ingredient in our construction is the function 

(fig^eic, z\h) = eg{e2C26)(l){{eiCi + Te2C2)9 + h,z), h e Er . 

Note that, in contrast with the classical r-matrix (4.11), there is an additional term 
ipg,,{O,Z-w\h)To0To. 

Proposition 6.1 The matrix R satisfies the quantum YB equation 

Ruiz - w)Ru{z)R23(.w) = R23{w)Rl3{z)Rl2{z - w) . (6.1) 
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Proof. 

Consider the coefficients in front of Ta'^T^'Si in the l.h.s of (6.1). It vanishes if a + 6 + c 7^ 0. 
Therefore (6.1) imphes 

ee{a X 6 + 26 X c)(pe,e{c, z - w\h)ipe,e{a - c, z\h)ipe,e{b + c, w\h) 

-ee{a X 6 + 26 X f)ipe,e{f, z - w\h)(pe,e{a - f, z\h)(pe,eib + f, w\h) = 

Here a, b are fixed and c, / are arbitrary elements of the lattice ^^^^(r) (C.16). Choose / in the 
form / = a — 6 — cto obtain 

ee(a X 6 + 26 X c) {<pff^^{c, z — w\h)(fo^^{a — c, z\h)(fo^^{b + c, w\h) 

-(fe^eia -b - c,z - w\h)(pg^^{a - c,w\h)(pg^eib + c, z\h)) = 0. 
This equality can be transform to the form 

ipe,e{cL + 2^, z\h)ipe,e{b, w\ti)ee{a x 6 + 26 x c) 

X (0,.(c + fi)- Ce,e{a -b-c + h)+ Ce,,(a -c + h)- Ce,,(6 + c + h))=0 

by means of (A. 27). If 6 = then the l.h.s. vanishes and we come to (6.1). Let 6 7^ and 
consider the shift c ^ c + jb, where j G Z. The shift does not change the exponential factor. 
Take the sum over the orbit generated by the shifts 

Co Ac + jb + h)- Ce,e(a -b-jb-c + h) + Ce,e(o -c + jb + h)- Ce,e{b + c + jb + h). 

One can see that the neighboring terms in the series vanish and we come to (6.1) for an arbitrary 
6. □ 

This Proposition allows us to define the associative algebra Pe,e,T,h by the relation 

R{z - w)L^{z)L^{w) = L^{w)L^l{z)R{z - w) , (6.2) 

where 

Lf'iz) = So^eeiO, z\h)To + ^ S.a^eAo^, z\h)Ta . 

a 

In order to obtain the relations in Pg ^-rh one should exclude the spectral parameters z, w from 
(6.2). 

Lemma 6.1 The relations in the associative algebra Pe,T,h assume the form 

Sa-cSb+c^e{c ■ (g - 6))f|9,e(a, b,c\h) = , for any a, 6 G Z © Z , (6.3) 

c 

where 

feM b, c\h) = Ce,e(c + h)- Ce^a -b-c + h)- CeAb + c + h) + Cg^a -c + h) (6.4) 

^ ^'me,e{c^ + 2h)^eAP)^eAP - a + 27) 

'?e,6(a - 7 + h)i^eAf^ + 7 + fi)^dAP - « + 7 + fi^^eA-y + ^) ' 

and T^o^e, CoA-) ^'^^ constants (C.17), (C.18). 
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Proof. 

Consider in (6.2) the matrix element ?&. We come to the relation 

X Sa-cSb+c{^d,e{(^^ ^ - w\K)^e,e{a - C, z\h)ipe,e{b + C, w\h) (6.5) 

c 

-(pe,e{a -b-c,z- w\h)ipo,e{a - c, w\h)cp0{b + c, z\h)}ee{c • (a - 6)) = . 
The expression in the brackets { } is equal to 

ee{{a + 2h)z + {b + h)w){4>0^e{c, z — w\h)4>0^e{a - c, z\h)4>0^e{b + c,z- w\h) 

-~4'0,€{O' — b — c, z — w\h)(j)0^^{a — c,w\h)(pQ^^{b + c, z\h)} , 

where (l)g^e{c, z\h) = ^((eici + Te2C2)0 + h,z). Finally, by using (A. 27) we come from (6.5) to 
(6.4). □ ' 

We rewrite (6.3) in the form most resembling the original Sklyanin relations [14]. Let 2t(a, (3) 
(respectively, ©(a,/3)) be the operator of antisymmetrization (respectively, symmetrization) 
with respect to the permutations of two indices (a,/3). 



Proposition 6.2 The relations (6.4) are equivalent to the commutator relations 

^ ieA^,b,l\n) 



iSoJ.]--ES-A.Mb-.)^^^^ (6.6) 



[Sa, Sp] = 2l(a, (3) I [Sa+pSoeei-P ■ a) - SoSa+peg,,{P ■ a)) f'' ^' (6-7) 



+ XI {Sa-ySp+^eei-f ■ {a - (3)) - Sp+^Sa-ySgi-j ■ (a - /?))) 

7y^a,-/3 

where CQ{a,P) are sl(A?^, C) structure constants (B.8), and 



tg,e{a,l3,0\h) 



6{a,P) \ (Sa+pSoegi-P ■ a) - SoSa+peg^P ■ a)) ^' "|^.^ 

+ X (Sa-^Sp^^eei'y ■ (a - /?)) - Sp+^Sa-^ee{-i ■ {a - /?))) ^'gj^j 



> =0. 



Proof. 

The relation (6.3) can be rewritten in the form 

X (^Sa-cSb+cegic ■ (a - 6)) - Sh+cSa-ceg{c ■ (a - 6))j ig^a, b, c\h) = . 

c 

due to the equality f0^e(a, 6, c|/i) = —igAa,b,a — b — c\h). For the case a = this relation assumes 
the form (6.6). To come to (6.7) we single out the two terms with c = and put them in the 
left hand side. In the right hand side we first write down the terms with c = a and c = —b and 
obtain the required relations. □ 
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It would be interesting to describe a quantum version of the whole bihamiltonian structure 
for elliptic rotators on sin* or Ham(T'^). Posing the problem more generally, the bihamilto- 
nian recursion precedure of generating the conserved quantities from a Casimir function for a 
linear family of Poisson structures (see [19]) might have a quantum analog as an expansion 
of a central element for a linear family of associative algebras. The latter seems to be a very 
strong requirement on a pair of associative algebras and it would be very interesting to find any 
non-trivial example of this kind. Namely, one is looking for a pair of associative algebras, such 
that their mixed associator satisfies some consistency condition. This would provide the most 
straightforward quantization for a system with a bihamiltonian structure. 
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7 Appendix 

7.1 Appendix A. Elliptic functions. 

Here we summarize the main formulae for elliptic functions. Consider an elliptic curve 

Er = C/(Z + rZ) , q = e(r) = exp(27riT) . (A.l) 
The basic element in our consideration is the theta function: 

nez 



-2inz\ 

± - q jyL - q c - q c ) 

n=l 



1. The Weierstrass functions: 



a{z; r) = exp(r?z )^^7^^ , (A.3) 
where 

C(z;T) = 5,logi?(z;r) + 2r?(T)z, C(^; r) - i + 0(^3) . (A.5) 

P{z;t) = -d,C{z;T) . (A.6) 

= ^ + + E ( + fcr + - (JtW) • ^^-^^ 
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2. Function 6: 



It has a pole at z = and the expansion 

cPiu, z) = - + C(«|r) + 2r]{T)z + |((C(tx; r) + 2r]{T)zf - p{u)) + ..., (A.9) 
as z ^ 0. Furthermore, 

4>{u, z)-'du(j)iu, z) = C(u + z)- C(u) + 2r?(r)z . (A.IO) 

(/>(n, z) = exp(-2ry«z)-^|^^-i-^ , (A.ll) 

^{u, z)^{-u, z) = p{z) - p{u) . (A.12) 

3. Quasi-periodicity: 

T?(z + l;r) = -??(z;r), + r; r) = -e(-^r - r) . (A.13) 

C(z + l;r) = C(^;r)-27?, C(z + r; r) = C(^; r) - 2(7ri + tjt) . (A.14) 

p(u + l;r) = p(« + r;r) = p(?x;r) . (A.15) 

(l){u + 1, z) = (f){u, z) , (f){u + T, z) = e{—z)(f){u, z) . (A.16) 



4. Parity: 



^(-2;r) = -i?(z;r). (A.17) 

C{-z-t) = -C{z;t). (A.18) 

P(-«;t) = p(n;T). (A.19) 

(f>{-u, -z) = -(f>{u, z) . (A.20) 



5. The Fay three-section formula: 



(f){ui, Zl)(p{u2, Z2) - (p{ui + U2, Zl)(p{u2,Z2 - Zl) - (p{ui + U2, Z2)(p{ui, Zl - Z2) = ■ (A.21) 

Particular cases of this formula are (A.12), the Calogero functional equation 

(j){u, z)dy(l){v, z) - (j){v, z)du(l){u, z) = {p{v) - p{u))(j){u + v,z), (A.22) 

and 

(t){ui,z)(t){u2, z) - ^(ui + U2,z){C,{ui) + C(u2) - 277(T)(ni + U2)) + dz^{ui + U2,z) = Q. (A.23) 
For ui + U2 + = 0, 

^{UI,Z)^{U2, Z)(I){U3,Z) = [p{z) - P{U3)] [({ui) + C{u2) + C(^^3 - z) + C{z)] , (A.24) 
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as follows from (A. 10), (A. 12), and (A.23). Then 



(I){UI,Z)<I){U2, Z)4){U3, z)\z^o = ^ + ^ [C(^l) + CM + C(^3)] 



(A.25) 



— p'M - Pius) [C(ni) + C(n2) + C{u3)] + 0{z) . 
^From (A. 21) and (A. 10) we have 

z)(l){u - ui, z) = (f){u, z){C{ui) + C{u - ui) - C{u + z) + C{z)) . (A.26) 
Another important relation is 

(p{v, z — w)(p{ui — V, z)(p{u2 + v,w) — (p{ui — U2 — v,z — w)(p{u2 + V, z)(p{ui — v,w) = (A. 27) 

(t){ui,z)(t){u2,w)f{ui,U2, V) , 



where 

f{ui,U2,v) = C{v) - C{U1 -U2-V) + ({Ul -V)- ({U2 + v) . 

One can rewrite the last function as 

i)'{0)i){ui)'&{u2)'&{u2 -ui + 2v) 



(A.28) 



(A.29) 



"^^ ^' ^' ^ 'd{ui — v)'d{u2 + v)'d{u2 — ui + v)'d{v) 
To prove (A. 27) one should consider the expression 

(f){v, Z — w)(f){ui — V, Z)(I){U2 + V,w) — (f){ui — U2 — V,Z — w)(p{u2 + V, z)(l){ui — V, w) 

(l){ui,z)(f){U2,w) 

It is a doubly periodic function in z and w without poles. Therefore, it is a constant, which 
depends on ui,U2,v. This constant is equal to f{ui,U2,v) (A.28). It is easy to check that the 
elliptic functions (A.28) and (A.29) coincide. 

A particular case of (A. 27), which corresponds to the case v = ui (or v = —U2), is the Fay 
identity (A. 21). Another particular case comes from ^2 = (or ui = 0): 



(l){v, z — w)(f){u — V, z)(f){v, w) — (f){u — v,z — w)(f){v, z)(j){u — V, w) 
^{ui,z){p{v) -p{u-v)). 

7.2 Appendix B. Elliptic constants related to gl(iV, C). 

Let 1 < k < N he & coprime number with respect to N and set 9 = j^. Define 

eg{z) = exp{2TTi6z) , 

Q = diag(e0(l), . . .,eg{m), . . . , 1) 
/ 1 ••• \ 



(A.30) 



A 



1 




V 1 







1 



(B.l) 
(B.2) 

(B.3) 
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Consider a finite two-dimensional lattice of order A'"^ 

1%^ ='L/N'L®'LIN'L. 

The matrices 

Ta = ^ee(^)g»iA»^ , (a = (ai,a2) G Z(2) (B.4) 
generate the basis in gl(A'^, C). We use the Greek letters for the elements of the lattice 

Zg)=zg)\(0,0). (B.5) 
Thus, {Ta} define a basis in sl(A/", C). Since 

Tan = ^ee{-'^)Ta+b, iaxb = aib2-a2bi) (B.6) 
the commutation relations in this basis assume the form 

[Ta,Tp] = Ce{a,p)Ta+p, (B.7) 

where 

C0{a,P) = —simreiax (3). (B.8) 

Let 

zfHr) = {-n + -f2T)e . (B.9) 

be a regular lattice of order iV^ — 1 on Introduce the following constants on 'L^q\t): 

Ml)=^{{li+l2r)e), (B.IO) 

Ce(7) =C((7i + 72t)0), Peil) = Piini + l2T)e) , (B.ll) 
and the quasi-periodic functions on 

Ml,z)=cl){iji+j2T)9,z), (B.12) 

veil, z) = 6^(722)^0(7, z) , (B.13) 

fed, z) = ee{-i2z)du(t){u, 2;)|„=(^j+^2^)0 . (B.14) 

It follows from (A.8) that 

Ve{l , z + I) = ee{'^2)ve{l , z) , (pe{i , z + T)eQ{--ii)ip0{'-i , z) . (B.15) 

7.3 Appendix C. Non-commutative torus 

1. DeGnition and representation. 

The non-commutative torus Ae is a unital algebra with the two generators {Ui, U2) that satisfy 
the relation 

U1U2 = e^^U2Ui, ee = e^-^^ 9 G [0, 1) . (C.l) 
Elements of Ae are the double sums 

A0=lx= Ca„a,U^'U^^ I Ca„a, € C i , 
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where Ca^_a2 elements of the ring © of rapidly decreasing sequences on 1? 

6 = {cai,a2 I sup„j^„2g2(l + «? + «i)*'|cai,a2p < ^ forall G N} . (C.2) 
The trace functional tr(x) on Ae is defined by 

tr(X) = coo • (C.3) 

It satisfies the evident identities 

tr(l) = 1, tr(XF) = tr(rX) . 
The dual space to 6 is the space 

01,02 

The relation with the commutative algebra of smooth functions on the two-dimensional torus 

= {M^/z ^ {0 < xi < 1, < a;2 < 1}. (C.5) 
comes from the identification 

C/i^e(x), [/2^e(y), (C.6) 
while the multiplication on becomes the Moyal multiplication: 

(/ * := /5 + E ^^n.i • • • (C.7) 

n=l 

The trace functional (C.3) in the Moyal identification is the integral 

tr/ = ^ fdxidx2 = foQ ■ (C.8) 

We can identify Ui, U2 with matrices from GL(oo). Define GL(cx)) as the associative algebra 
of infinite matrices CjkEjk, where Ejk = \\Sjk\\, such that 

suPj-fegz|cifePli - fcl" < 00 forn G N. 

Consider the following two matrices from GL(oo): 

Q = diag(e(j6')) and A = , j G Z. 

We have the following identification 

Ui^Q, U2^A. (C.9) 
Another useful realization of Ae in the Schwartz space on M by the operators 

Uif{x) = f{x-e), i72/(a:) = exp(27rix)/(x) . (C.IO) 

2. .sin-algebra. 

Define the following quadratic combinations of the generators 

r. = ^e(^^)TO aGZ(^ (C.ll) 
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= {a = {ai,a2),aj G Z, a 7^ (0,0)}. 
Their commutator has the form of the sin-algebra 

[Ta,Tfs] = Ce{a,P)Ta+p, (C.12) 

where 

Ce{a,(3) = ^simreiax P) . (C.13) 
We denote by sing the Lie algebra with the generators over the ring © (C.2) 

^ = J2i'aTa, i^aee, (C.14) 
a 

and by SINg the group of invertible elements from Ao- The coalgebra sm| is the linear space 



S = ^ SaTa, Sa e 6' 



sing = < 



In the Moyal representation (C.7) the commutator of sing has the form 

[f{x, y), g{x, y)] = {/, g}' := lif * g - 9 * f) (C.15) 



3. Elliptic constants related to NCT Ag 

Introduce two numbers e = (ei, €2) such that < 1 and eaO are irrational. Consider the dense 
set Z0,e(r) in Et-: 

Ze,6(r) = {(ei7i + re272)^e^r I (71, 72) e Z^^)} . (C.16) 
The corresponding elliptic functions with the arguments from Z5t_e(r) are as follows: 

%(7) = ^?((ei7i + e272r)e) , (C.17) 

Ce,e(7) = C((ei7i + ■^€272)6') , pg,eil) = piieiji + Te2j2)0) . (C.18) 

<^0,,(7, z) = <^(-(ei7i + re272)^, z) , (C.19) 

feAl) = M^2l2z)(i)g,e{7, z) . (C.20) 

fe,e{l) = eg{€2j2z)du(t>{u, Z)|„=(ei7i+r6272)e " (^-21) 

4. Dispersionless limit. 

In the limit ^ — > the Lie algebra sing becomes the Lie algebra of Hamiltonian functions 

Ham(T^) ~ C°°{T^)/C (C.22) 
equipped with the canonical Poisson brackets. In Ham{T'^) we have the Fourier basis 

e(a • x) = exp(27ri(aixi + 02X2)) (C.23) 
instead of the basis (C.ll). The commutator (C.12) becomes 

[e{ax),e{l3x)] = (a x P)e{{a + P) ■ x) . 
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The algebra Ham{T'^) (without constant Hamihonians) is isomorphic to Lie algebra SVectQ(T'^) 
of the divergence-free zero- flux vector fields on equipped with the area form —A7r^dxidx2- 
Let h{xi,X2) e Ham{T^). Then the Hamiltonian field Vh corresponding to the Hamiltonian 
function h is 

Vh = -^iid2h)di - {dih)d2) , (C.24) 

while 

[^^,^/.'] = (C.25) 

For f{x) = fae{a ■ x) 



,.-f = -S^A. (C.26) 
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